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Origin of "Very-Low-Frequency Emissions" ' 

R. M. Gallet and R. A. Helliwell 2 

(January 16, 1959) 

Selective traveling-wave amplification in the outer ionosphere is postulated to explain 
very-low-frequency emissions, a class of very low-frequency (1 to 30 kilocycles per second) 
natural noise. By analogy with the mechanism of traveling wave tubes, low-level ambient 
noise in the outer ionosphere is amplified in streams of incoming ionized solar particles at 
frequencies for which the stream and wave velocities are equal. Required velocities are 
in the range 0.01 to 0.1 c (where c is the velocity of light). Streams with densities of the 
order of one electron per cubic centimeter would provide sufficient energy. Phenomena 
which can be explained qualitatively by the theory are the hiss, quasi-constant tones, dawn 
chorus and related transients, and very long trains of whistler echoes. A quantitative 
example shows how the theory can reproduce the general form of certain characteristic 
discrete spectra "hooks" of emissions, and how this leads to definite values of particle 
velocity and a law for the distribution of electron density in the outer ionosphere. 



1. Introduction 

When a high-gain audio amplifier is connected 
to a loop or long-wire antenna, unusual sounds 
of natural origin are often heard. The more 
impulsive sounds called "atmospherics" are pro- 
duced by lightning and reach the receiver by 
traveling between the ground and the lower edge 
of the ionosphere. They are described by various 
terms such as "cricks," "tweeks," and "crashes." 

Some of the energy from a lightning discharge 
also enters the ionosphere and is guided by the 
lines of force of the earth's magnetic field. As 
the energy propagates from one hemisphere to the 
other, it is dispersed and arrives at the receiver 
as a whistling tone, usually of descending frequency. 
Such signals are called "whistlers," and their 
theory has been developed in some detail (Storey 
[1] , 3 Helliwell et al. [2], Ellis [3], and Maeda [4]). 

There is yet another class of audiofrequency signals 
which we will call "very-low-frequency emissions". 
These appear to be related to magnetic disturbance 
and occur most frequently near the auroral zones. 
(Storey [1], Pope [5], Allcock [6], Dinger [7], and 
Watts [8].) One of the most common is the "hiss", 
a band of noise several kilocycles wide which is 
usually heard below 10 kc; it is occasionally ac- 
compained by one or more discrete tones in the same 
range of frequencies. Another is the so-called 
"dawn chorus", a combination of short (0.1-0.2 
sec) rising whistles and warbling tones. Various 
isolated events, including rising and falling tones and 
combinations thereof can be observed. Very-low- 
frequency emissions are frequently associated with 
strong whistlers. Sometimes a whistler appears 
to " trigger" short periods of dawn chorus or hiss. 
Whistlers are sometimes preceded or followed by 
one or more rising tones. 

It is the purpose of this paper to outline a theory 
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of the origin of vlf emissions, based on selective 
traveling-wave amplification taking place in the 
outer ionosphere. The input signal is assumed to 
be provided by thermal radiation, whistler energy, 
or possibly Cerenkov radiation (Ellis [9]). It is 
assumed that the energy for amplification of the 
electromagnetic wave is provided by the streams 
of ionized solar particles (both positive and negative) 
which are thought to cause the auroras. These 
streams are assumed to penetrate the ambient 
ionization of the outer ionosphere and to be guided 
along the lines of force of the earth's magnetic field. 
The mechanism of amplification is essentially similar 
to that which takes place in an ordinary traveling 
wave tube ( Tierce [10]), except that the "slow-wave" 
circuit (provided by the helix in the tube) is the 
ambient ionization of the outer ionosphere in the 
presence of the earth's magnetic field. Such a 
medium is dispersive, i.e., the wave velocity is a 
function of frequency, and at the vcvx low fre- 
quencies considered here the velocity may be reduced 
as much as two orders of magnitude below that in 
free space. Furthermore, as shown by Storey [1] 
in connection with the theory of whistlers the energy 
of a vlf wave tends to be guided along the earth's 
magnetic field lines rather closely, and thus may 
follow the same path as the postulated particle 
stream. The wave normal, on the other hand, may 
make any angle with the field; and since the wave 
polarization is approximately circular, there will 
normally be a substantial component of the electric 
field of the wave which is parallel to the magnetic 
field. This component will act on the moving particle 
stream to produce the required sinusoidal modulation 
of charge density. For given values of the ambient 
plasma and gyro frequencies, the stream velocity 
will approximate the phase velocity at two par- 
ticular frequencies, and traveling wave amplification 
can take place. As the parameters of the medium 
change with position along the stream, the fre- 
quencies of amplification will change. This paper 
develops the conditions for interaction assuming 
the validity of the analogy with the traveling wave 
tube mechanism. 
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2. Interaction Conditions 

By analogy with the small-signal theory of 
traveling wave tubes (Pierce [10]), the necessary 
condition for coupling between the electromagnetic 
wave and the stream of charge is 



v v =v S: 



(1) 



where v p = phase velocity of wave, v s = velocity of 
stream. 

We must now deduce the conditions under which 
(1) will be satisfied. The phase velocity is en, 
where c is the velocity of light and n is the refractive 
index in the medium. Since we need a phase velocity 
less than the velocity of light, we shall be interested 
in the "whistler-mode" of propagation (Storey [1]). 
For simplicity we shall consider the case in which 
the wave normal direction is approximately the 
same as that of the earth's magnetic field. The 
refractive index for the propagating mode is then 
given by 



n 6 - 



= !+/!/{/(/*-/): 



and the phase velocity is then simply 
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(2) 



where j=w&re frequency; /#=gyrofrequency; /,,= 
(8 IN)M= plasma frequency; with N= number of 
electrons per cubic centimeter when/,, is in kilocycles 
per second. 

It is assumed that /„>/,/»>/. From (1) and (2), 
the frequencies of amplification are given by 



1± 
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If we restrict the discussion to cases in which the 
stream velocity is less than about 0.1 c, then (v s /c) 2 <^l, 
and (3) can be written 



/=^{i±[i-(2P//„)r 2 }, 



(4) 



where P—J p (v s /c) and is conveniently measured in 
kilocycles per second. 

For interaction to occur, the solution to (4) must 
be real, and hence the condition for interaction is 



f H ^2P. 



(5) 



The maximum value of the wave frequency / is j H , 
since the maximum real value of the radical in (4) 
is unity. The conditions for interaction are shown 
in figure 1, which is a plot of wave frequency/ as a 
function of gyrofrequency f H for different values of 
the parameter P. All curves are asymptotic to the 
line j=j H (ignoring the dashed curve for the mo- 
ment). 



As an example of the application of these curves, 
consider a stream of velocity 3,000 kmps which enters 
the outer ionosphere along an auroral line of force 
(near geomagnetic latitude 67°). The ambient 
plasma frequency is assumed to be constant and 
equal to 200 kc, and thus the parameter P=3 kc. 
When the stream reaches the point where j H =2P, 
interaction can begin at the wave frequency /^/2, or 
3 kc. Beyond this point two frequencies can be 
amplified, and as j H increases the lower one ap- 
proaches zero and the upper one approaches f H , 
following the curve P=3 in figure 1. The locations 
of these frequencies along the dipole line of force 
are plotted in figure 2. 




Figure 1. Variation of excited wave frequency with gyro- 
frequency. 
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Figure 2. Illustration of the theory of hiss, for the auroral line 
of force X= 67°. 
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3. Stream Velocities 

The stream velocity required for 
found from (4) and is given by 



J V 



interaction is 



(6) 



Values for the stream velocity for appropriate 
values of the medium parameters and the wave 
frequency are given in table 1. It is seen that the 
required stream velocities could vary roughly be- 
tween 0.01 and 0.1 of the velocity of light. These 
velocities seem reasonable, since the lower value, 
0.01 c, is comparable to the lower limit found by 
Meinel [11] in spectroscopic studies of hydrogen in 
the aurora. On the other hand, it appears necessary 
that electrons which penetrate to the E region must 
have velocities of approximately 0.3 c (Vegard [12]). 
Thus the velocities required for the traveling wave 
mechanism appear to lie in the same range as those 
ascribed to auroral particles. 

4. Stream Energy Requirements 

If traveling wave amplification is to be accepted 
as a possible mechanism, it is necessary that the 
power density of the incoming stream exceed that 
of the observed signals. A measurement of the 
strength of relatively strong hiss was made at 
Stanford, California (Stanford University [13]) on 
March 3, 1956, during a period of magnetic dis- 
turbance. The power density was of the order of 
10~ 10 w/m 2 . for a signal having a bandwidth of about 
1 kc. Observations of whistlers indicate that the 
total attenuation from the source to the observer is 
not likely to exceed 20 db. Thus the power density 
of the signal as it leaves the amplifiying region is of 
the order of 10~ 8 w/m 2 . 

Let us now calculate the power density of the 
incoming stream. Since the particle densities and 
velocities of streams from the sun are not known, 
we must make a conservative estimate. From the 
work of Martyn [14] the stream is estimated to have 
a density of at least 1 electron/cm 3 and a velocity of 
3,000 km/s. The power density in such a stream is 
(%mv 2 )(Nv), where m is the particle mass, v the 
velocity, and N the number of particles per unit 
volume. For electrons, m=9.11X10~ 31 kg, and 
the power density of the stream is thus about 



10~ 5 w/m 2 . The required efficiency of conversion 
from kinetic energy to electromagnetic energy is 
thus 0.1 percent, a low value compared with traveling 
wave tube efficiencies, which are of order 10 percent. 

5. Comparison of Theory With Observations 

Using the traveling wave hypothesis the hiss is 
explained in terms of a continuous stream of particles 
along which amplification occurs only for a restricted 
band of frequencies. The amount of gain at any 
one frequency will depend upon the interaction 
distance available, which in turn will depend on the 
way in which the plasma frequency varies with dis- 
tance along the path. In a like manner the quasi- 
constant tones sometimes observed on different fre- 
quencies can be explained by assuming an ionization 
model in which the plasma frequency increases with 
gyrofrequency in such a way that /= constant over 
an appreciable length of path. Several different tones 
observed simultaneously could be explained by 
different streams flowing along lines of force of 
sightly different geomagnetic latitude, or possibly 
by streams of different velocities. 

Transient phenomena, of which the dawn chorus 
is an example, are less easily explained. The travel- 
ing wave hypothesis can be adapted to the generation 
of transient signals by assuming that the incoming 
particles arrive in bunches. As the bunches travel 
along the lines of force, amplification takes place in 
the region of the bunch. However, the character- 
istic frequencies are amplified for a length of time 
which depends upon the length and velocity of the 
bunch. As the bunch moves, the frequency which 
is generated changes in accordance with changes in 
the parameters of the medium. The mechanism of 
formation of such bunches is not known; however, 
the visual phenomena observed during aurora sug- 
gest that solar streams may be bunched. 

The interpretation of observed transient signals 
in terms of bunches of particles is complicated by 
the fact that the frequency-time relation of the 
received signal will depend on the integrated group 
delay from the point of origin to the receiver. The 
group delay, which is a function of frequency, 
clearly will affect the frequency- time dependence of 
the received signal. However, it cannot by itself 
explain the appearance of a given frequency at two 
different times, a phenomenon seen frequently in 
the dawn chorus. Such cases might better be ex- 
plained in terms of variations in the ambient ioniza- 



Table 1. Calculated stream velocities for different parameters of amplifying region 



Gyrofrequency 


Distance from surface 


Minimum geomagnetic latitude 


Plasma 
frequency 


Wave frequency 




2kc 


5kc 


10 kc 


Sh = 1,000 kc 


Regular layers 


Middle latitudes 


4, 000 kc 

1, 000 kc 

1, 000 kc 

200 kc 

200 kc 

200 kc 


0. 011 c 
.044 c 
.014 c 
.070 c 
.02 c 
.014 c 


0. 018 c 
.071c 
.022 c 
.11 c 
.025 c 


0. 025 c 


/// = 100kc 


'--'1.0 earth radii 


44° 


,1c 
.030 c 


/„=10kc-_ 


3.6 earth radii. _ -.. ._ 


61° 


.15 c 


0/=5kc __ 


9.0 earth radii .. _ _. 


65° 
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tion density in the interaction region. For example, 
consider the hypothetical model, / p =200 kc+/ H 2 / 
(1 kc), where j H is measured in kilocycles per second. 
A short bunch of charge (say, with dimensions of 
about 100 km) traveling through such a region 
would excite the frequencies given by the dashed 
curve of figure 1. The energy would be limited to 
short groups produced at different times. The 
curve of excited frequency versus gyrofrequency 
transforms to a somewhat similar shape in the 
frequency-time plane. Some distortion of this 
curve occurs at the receiver because group delay is 
a function of both frequency and location of source. 
No completely closed trace of this type has yet been 
observed, but many of the recorded signals first fall 
and then rise with time, corresponding to the lower 
portion of the dashed curve. Thus there is at least 
a qualitative resemblance between the theory and 
the experimental data. 

The steady hiss and the sharply defined transients 
represent extremes among vlf emission phenomena. 
In between lie many forms which show varying 
degrees of irregularity. Their character often 
changes markedly in the course of a few hours. It 
is suggested that the observed variety of form is 
simply the result of variations in the degree of bunch- 
ing of the charge in the streams. 

The traveling wave hypothesis offers a convenient 
explanation of very long trains of echoes of a 
whistler. The echo intensity decreases very slowly 
with order, and may even increase. It is difficult 
to explain this behavior in terms of the passive 
properties of the medium. Discrete echoes in a 
single train have been detected for as long as 4 min. 
These rather rare events are usually accompanied 
by magnetic disturbance and the dawn chorus or 
hiss. The whistler components which are most 
prolonged often lie in the 3- to 4-kc range where 
the hiss is known to predominate. On the basis of 
these facts it is suggested that the low apparent 
attenuation in very long whistler trains can be 
attributed to traveling wave amplification some- 
where along the path. The whistler energy propa- 
gates along the same path followed by the incoming 
solar stream and is amplified in the regions where the 
whistler phase velocity approximates the stream 
velocity. 

One of the most puzzling phenomena is the inter- 
action between loud whistlers and vlf emissions. 
If during a period of dawn chorus and whistler 
activity, an ordinary loud whistler should occur, it 
is frequently accompanied by one or more rising 
whistles (Storey [1]). In attempting to explain this 
connection in terms of the traveling wave concept, 
we suggest that the loud whistler may be capable of 
accelerating an incoming stream of particles. For 
energy to be transferred from the whistler to the 
stream, the wave must accelerate the particles in 
its electrical field. Equation (2) shows that in a 
region where j v does not increase too rapidly with 
j H , the phase velocity of a fixed frequency wave will 
increase with j H . For our purposes the rate of 
change of whistler frequency with time is sufficiently 



small that we can assume the accelerating signal to 
be of constant frequency. Thus as the whistler 
travels down the field line, particles of a given 
velocity are literally picked up by the wave at every 
point where a wave component of the same velocity 
exists. As the wave continues to propagate, its 
phase velocity increases and the particles riding on 
the front of the wave are accelerated, as in a linear 
accelerator. With the additional kinetic energy 
which the stream acquires from the signal it would 
be more efficient as a traveling wave amplifier and 
might conceivably produce audible vlf emissions. 

6. Conclusions 

The rather close similarity between conditions in 
the outer ionosphere and those in a traveling wave 
tube lead to the hypothesis that traveling wave 
amplification might take place at very low audiofre- 
quencies. Even though the traveling wave hy- 
pothesis may not account successfully for all of the 
observed phenomena, it is probably an important 
factor to be considered in whistler propagation as 
well as in the generation of vlf emissions. Confirma- 
tion of the hypothesis must await a more detailed 
study of the traveling wave-gain problem and some- 
what more detailed information regarding the dis- 
tribution of ionization in the outer ionosphere. Such 
data presumably will come from whistler analysis. 
This theory, if correct, should provide a useful new 
tool for a study of the densities, velocities, and loca- 
tions of streams of charged particles in the outer 
ionosphere. 

7. Appendix: A Quantitative Comparison 
of the Theory With Observation 

7.1. Example Chosen for Comparison With the 
Theory 

The theory presented above has been developed in 
extensive numerical applications by one of us 
(R. M. G.) to see if the theory is able to account in 
detail for spectra of discrete vlf emissions, using a 
minimum of reasonable assumptions. 

There is considerable information on the distribu- 
tion of electron density up to the height of maximum 
density in the ionosphere, but very little is known 
about densities above this height. In particular, the 
law governing the distribution of electron density 
with height is not known for this region. In these 
investigations several models of the distribution 
above the main body of the ionosphere were con- 
sidered. They were: Constant electron density; a 
power law r~ n where r is the distance from the center 
of the earth and n was taken from 1 through 6 ; and 
an exponential distribution e~ kT where k is a con- 
stant. The power law was found to give the best 
representation of certain distinguishing features of 
the vlf emissions and also to be consistent with in- 
formation obtained from parallel studies of whistlers 
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(dispersion constants and analysis of nose whistlers). 
The same criteria seem to require the exponent 
nc^'A. Furthermore, a power law distribution with 
nc^S appears to be justified theoretically in other 
research in progress on the gravitational equilibrium 
of extended atmospheres, at large distances from the 
earth's surface. It is perhaps significant for the 
physics of the exosphere that nc^3 menus that the 
electron density is nearly proportional to the strength 
of the earth's magnetic field. 

An example of a computed spectrum (frequency 
versus time at the earth's surface) of waves emitted 
by a discrete cloud of ionized particles moving with 
constant speed parallel to the lines of force of the 
earth's magnetic field is shown in the upper part of 
figure 3. This example was chosen to approximately 
match the observed spectrum shown in the lower 



part of the figure. The latter is a typical member 
of the "Hook" class of vlf emissions. This class is 
one of the best defined, and the reproducibility of 
the emissions is remarkably good in general shape. 
However, no two hooks, even though separated by 
only a few seconds in time, are exactly the same. 
Therefore, the exact matching of any individual tone 
is of no particular importance. In the example 
shown here, no great accuracy in matching was at- 
tempted; the same computations made with a 
slightly lower speed for the cloud of particles (say 
9,000 km/sec) would have resulted in a better match. 
It is more significant to realize that with reason- 
able numerical values for the parameters, the theory 
leads very naturally to calculated spectra very 
similar to those observed in regard to values of 
frequency, duration, and general shape. 
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Figuke 3. Comparison of theory and observation {see appendix for details). 
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7.2. Model of the Distribution of Electron Density 

The model of electron density versus altitude used 
for the numerical calculations is illustrated in figure 
4. Electron density (iV e /cm 3 ) is given on the right 
scale, and the corresponding local plasma fre- 
quency i p in kilocycles can be read on the left scale. 

In the exosphere, far above the earth's surface 
(above 1,000 km), rather than r~ 3 , the distribution 
chosen was: 



NelH= constant, or f 2 p =Kf Hj 



(7) 



where H is the local magnetic field strength and f H 
is the local electron gyrofrequency. Below 1 ,000 km, 
as a means of joining smoothly with a plausible 
ionospheric distribution, the above expression for 
j% is multiplied by the factor exp{ — (h— 1,000) 2 / 
Jfex}, where A is the height in kilometers above the 
earth's surface. The scale height J^ x is taken as 
200 km. 

The ionospheric distribution is taken to be a 
Chapman-like layer 

JVa=iV max exp {^(1-z-O, 
or 

Pp,iono=f 2 p,max exp {^(l — Z — e" 1 )}, ( 8 ) 

where z= (h—h m&x )/j^ The chosen value of 2V max (10 6 
cm" 3 ) corresponds closely to the critical frequency 
of the F2 layer at the place and time of the observa- 
tion to be matched (Boulder, Colorado, October 26, 
1956, 2230 UT or 1530 local standard time). The 
value of h m&x was taken as 350 km. Above h m&x the 




1000 
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2000 
EARTH'S SURFACE, 



Figure 4. Distribution of electron density used for the theoretical 
calculation shown in figure 8. 



value of the scale height used was j^ = 200 km. 
Below h max , j^=100 km was used. 

Thus the total electron density for the model used is 

/2 {2 I /2 

J v J pjlonoTJ pjexo 

=/!,max exp {mi-z-e")} 



xfKSn exp {-(A-l.OOO) 2 /^} for A<1,000 km\ 
\Kf H for h> 1,000 km J 



+ 



(9) 



Note that only a small range of values of the con- 
stant K is possible if the electronic density in the 
exosphere is not to be too large near the earth (smaller 
than the lowest values of iV max ) nor too small to per- 
mit the propagation of whistlers at the greater 
distances. The value chosen here was 4,000 kc. 

7.3. Geometric Conditions 

The calculations were made for an ionized cloud 
of particles moving along the line of force which in- 
tersects the earth's surface at the geomagnetic lati- 
tude A=55°. The observations were made at 
Boulder at geomagnetic latitude 49°, but it is known 
from studies of whistlers that a signal propagating 
in the whistler mode can be received at places far 
from the foot of the nominal line of force along which 
the wave is propagated. For a dipole field, the line 
corresponding to X=55° intersects the plane of the 
geomagnetic equator at 2.04 earth radii, or 13,010 km 
above the earth's surface. (Note that fig. 4 shows 
the model only out to the height 3,500 km.) The 
total length along the line of force from the equator 
to the earth's surface is 3.17 earth radii, or 20,214 
km. At the equator the magnetic field strength 
along the line of force is a minimum. Its value is 
0.011 gauss compared to 0.315 gauss at the surface. 
The corresponding gyrofrequency is 31.4 kc, and 
with the exosphere model (7) N e = 1,560 and/ p =354 
kc at this height. 

Using as the origin of time a starting position in 
the equatorial plane, the calculations were made 
for points regularly spaced along the trajectory to 
obtain: (1) the frequencies emitted; (2) the time 
of the emission; (3) the time for propagation of the 
wave from this point to the earth's surface along 
the remaining arc, using the full formula for group 
velocity in the whistler mode; and, (4) the sum of 
these two times to obtain the arrival time at the 
earth. 

At 10,000 km/sec the total time of flight for the 
cloud is nearly 2 sec while the range of arrival times 
of the emitted waves is about % sec. Of course, 
one observes only the range of arrival times. 

The "knee" in the "hook" in the present example 
corresponds to the transition at about 2,500-km 
height between the electron density characteristic 
of the exosphere and the more rapidly increasing 
density in the ionosphere, which is responsible for 
the high branch of the "hook." 
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7.4. Velocity of the Cloud 

Generally speaking, the shape of the computed 
spectrum will remain approximately the same for 
different velocities. But the frequency and time 
scales will vary rapidly. A small velocity furnishes 
a range of arrival times that are too long and emitted 
frequencies that are too low. The durations and 
frequencies are very sensitive to the velocity. 
Velocities of about 1,000 kin/sec, as suggested by the 
delay between related solar and terrestrial events, 
give frequencies roughly ten times too small and 
durations 10 times too long. Satisfactory fits are 
obtained only for velocities of the order of 10,000 
km/sec. Such a high velocity seems to indicate that 
some process in the interaction between the earth's 
exosphere and the solar corpuscular radiations gives 
high local velocities to the small clouds responsible 
for the vlf emissions. These clouds are tentatively 
associated with the flashes of luminosity in active 
auroras. 

It is significant that this velocity is of the same 
order as that required for protons to penetrate to 
the 100-km level in the aurora. 

Table 2. Summary of the conditions assumed for the sample 
theoretical calculation 



Distribution of electron density 


Other conditions 


Ionosphere 


Exosphere 


X = 55° 

r, = 10< km/s 


iVm«=10»/cm« 

Am ax =350 km 

//, = 200km 

J^B = 100km 


Aex = 1,000 km 
jC = 200km 
K= 4,000 kc 
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